Are we seeing the beginnings of Inflation? 



On 
O 
O 
(N 

< 



O 

u 

Or 
i 

o 
in 

o3 



On 
Os 
O 

00 

o 

OS 

o 



- 1—^ 

x 



Cosmin Ilic 1 ,Tirthabir Biswas 2 and Katherine Freese 1 
1 Michigan Center for Theoretical Physics, Physics Dept., 
University of Michigan, Ann Arbor, MI 48109, USA 
2 Department of Physics, Pennsylvania State University, University Park, PA, 16802-6300, 

(Dated: August 7, 2009) 



USA 



Phantom Cosmology provides an unique opportunity to "connect" the phantom driven (low en- 
ergy meV scale) dark energy phase to the (high energy GUT scale) inflationary era. This is possible 
because the energy density increases in phantom cosmology. We present a concrete model where the 
energy density, but not the scale factor, cycles through phases of standard radiation/matter domi- 
nation followed by dark energy/inflationary phases, and the pattern repeating itself. An interesting 
feature of the model is that once we include interactions between the "phantom fluid" and ordinary 
matter, the Big rip singularity is avoided with the phantom phase naturally giving way to a near 
exponential inflationary expansion. 



PACS numbers: 98.80.Cq,11.25.-w, 04.50.+h 



I. INTRODUCTION 

The current accelerated expansion of the universe is 
usually explained by invoking a Dark Energy (DE) com- 
ponent 1 which today comprises more than 70% of the 
total energy in the universe (for reviews see 0, [H, [||). 
The case of a pure cosmological constant, with wa = 
Pa/ Pa = — 1 marks the divide to the 'phantom' realm. 
Phantom dark energy models are described by systems 
with 



&<-l 
Pp 



(1) 



and have the intriguing feature that the energy density 
in the universe increases with expansion, 



Pp 



,-3(1- 



(2) 



Hence a universe with low ~ mev scale accelerated ex- 
pansion can eventually reach energy scales close to the 
GUT scale, for instance. For some examples of cosmo- 
logical scenarios using phantom energy see [H, 0]. The 
question that we want to ask is whether it is possible to 
exploit this feature of phantom cosmology and turn the 
dark energy driven acceleration into a GUT scale infla- 
tionary phase 2 . The idea then would be to construct a 
cyclic model where dark-energy /inflationary phases are 
interspersed with decelerating radiation/matter phases. 

Several problems immediately appear. Firstly, unless 
the equation of state for the phantom phase, w p , is ex- 
tremely close to —1, the phantom acceleration will be 
much faster than the deSitter expansion, and cannot be 
reconciled with data. Density perturbations produced 



by 



1 For alternative approaches which try to avoid dark ener; 
invoking large scale inhomogeneities see, for instance, Qj, [2 

2 For a brief list of papers that propose various mechanism of con- 
necting the current accelerated expansion to inflation see 



during a phantom phase will give rise to a blue spec- 
trum, and consistency with the current WMAP 5-yr data 
at the 2a level with tensor modes included 0] requires 
— 1 > w p > —1.01. Secondly, it is well known that 
phantom cosmology typically ends in a Big-rip singu- 
larity, rather than the standard radiation phase which 
follows inflation. Remarkably, we find that both these 
problems can be addressed when we include interactions 
between the "phantom fluid" and some hidden sector 
matter. Such interactions ameliorate the phantom ac- 
celeration phase to an asymptotic de Sitter type expan- 
sion, once the phantom energy density reaches a criti- 
cal value. It is easy to arrange this transition to occur 
around the GUT scale, which is appropriate for inflation- 
ary cosmology. This also automatically avoids the big rip 
singularity as the space time now approaches a deSitter 
universe. The density perturbations can have a variety 
of possibilities, allowing for agreement with observations 
Q. Moreover, in our scenario the universe transitions to 
an asymptotic deSitter phase independent of the value of 
Up, and thus avoids having to fine-tune w p very close to 
— 1. As an additional advantage over the usual slow-roll 
inflationary scenario, in our phantom-driven inflationary 
model one does not have to tune the flatness of the po- 
tential usually necessary to obtain the large number of 
cfoldings and near scale-invariant spectrum. In addition, 
the hierarchy between the ~ meV dark energy scale and 
the GUT inflation scale can re-expressed in terms of four 
parameters that take values of 0(1) to 0(1O 2 ). We can- 
not however address the "coincidence" problem in our 
picture. Finally, there is the question of how to construct 
a theoretically self-consistent model of phantom energy. 
We will comment on this problem shortly. 

Before delving into the details of our specific realiza- 
tion of the "phantom cyclic model" , let us outline the ba- 
sic picture by considering just a simple two fluid model, 
phantom matter (p p ) + radiation (p r ). The cosmology 
we want to realize is the following: although the scale 
factor always increases monotonically with time, the en- 
ergy density "cycles" , at least approximately. Each cycle 
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is divided into two different phases: (a) Radiation domi- 
nated phase, which starts at an energy density p r = A 4 ax . 
As the universe expands radiation gets diluted, the Hub- 
ble parameter decreases and reaches a minimum when 
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From here on we enter (b) the phantom 



energy dominated phase. In realistic cosmology the ra- 
diation phase should give way to matter domination at 
energy densities ~ (10 eU) 4 , before giving way to phan- 
tom domination, but for simplicity we are going to ig- 
nore this slight complication. Thus for a typical scenario 
which would be consistent with dark energy and infla- 
tionary paradigm, A max ~ 10 15 GeV ~ 10~ 3 M p or the 
GUT scale, and A min ~ meV ~ 10~ 30 M p corresponding 
to the scale of current energy density. Now, in the phan- 
tom phase, as the universe expands the energy density 
increases, and so does the Hubble rate. Initially, depend- 
ing upon how negative the phantom equation of state 
parameter, u> p , is this increase in energy density can be 
quite fast. However, in our model we will see that once 
the energy density reaches close to a critical scale A max , 
which is determined by the interactions between the hid- 
den and ordinary matter sector, the energy density and 
the Hubble parameter asymptote to a constant giving rise 
to a near exponential expansion. This inflationary phase 
can end via the reheating mechanism described in section 
Mil IIVI after which we enter the radiation dominated era 
of the next cycle. 

A similar idea to the one in this paper has previously 
been presented by [H| who dubbed this model "the eter- 
nally expanding cyclic universe." They too (in their Sec- 
tion 4.3) suggested an alternating increasing/decreasing 
energy density. However the theory behind their model 
is quite different from ours, and consequently their pre- 
dictions for resultant density perturbations are different 
as well 3 . One major problem of all phantom type mod- 
els is the vacuum stability due to the null energy condi- 
tion violation. [l3j examine a consistent way to solve this 
problem, based on a deformation of the ghost condensate 
model of Arkani-Hamed et all (l7| . 

Although the interactions between phantom fluid and 
ordinary matter can lead to an inflationary spacctimc, 
we are still left with a graceful exit problem, or how to 
ensure that the universe enters the standard radiation 
dominated era. Depending upon the specific model differ- 
ent "reheating mechanisms" may be able to trigger such 
a transition. We focus on a model where the phantom 
fluid consists of a ghost like scalar field coupled to some 
hidden matter sector. Such a fluid closely resembles the 
interacting DE-DM models [l§ M EL S except that 
the scalar field instead of being an ordinary quintessence 
field, has negative kinetic energy like a ghost. Although 
field theory with ghosts is plagued with problems of uni- 



3 They always predict a negative, even if extremely small, tilt of 
the spectral index while our model will allow a variety of possi- 
bilities (as discussed in later sections of this paper). 



tarity/instability [Tg| , recent developments attemptin g to 
address these problems include progress in non-local 24 1 
and Lee- Wick (2f| higher derivative models; see also 13 j 
and [3. As we will see, in our model the transition from 
phantom to radiation phase and vice-versa is achieved 
partly by suitably choosing the interaction strength be- 
tween the scalar field and the hidden matter sector, and 
partly due to the presence of interactions between the 
hidden and ordinary matter sector. (There is no direct 
coupling between the ghost field and ordinary matter.) 

As emphasized before, the main reason why the cos- 
mology described above can replace the standard infla- 
tionary paradigm is because in the phantom phase the 
energy density increases even though the universe con- 
tinues to expand. Thus after the usual dilution in a radi- 
ation dominated phase, the phantom phase followed by 
reheating ensures that the universe again becomes hot 
and therefore can reproduce the successes of the Big Bang 
Model, such as Big Bang Nucleosynthesis and Cosmic Mi- 
crowave Background Radiation. There is another essen- 
tial similarity between our scenario and the inflationary 
paradigm. The essential reason why inflation solves the 
standard cosmological puzzles is because our observable 
universe (of radius ~ Hq 1 ) can originate from a very tiny 
region at the beginning of inflation. Something very sim- 
ilar happens in our model as well, the universe expands 
by a huge factor in every cycle. In our model the number 
of e-foldings in the radiation and the phantom phase is 
given by 

7V rad ~ In ( ) an d 7V" phan ~ — - 4 — - In ( 

' (3) 

In order to have a successful GUT scale inflationary 
paradigm we need J\f± n f > 60. Thus one gets: 
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Just to get an idea, if we take A m i n 

10 15 Gew, uj p 1.3, and Af ini ~ 60, we get Af tot ~ 400. 

What this means is that only a very very tiny portion (for 
the chosen example, e -400 th) of our observable universe 
will ultimately grow to become the observable universe in 
the next cycle at the same energy density scale. Another 
essential similarity between the standard inflationary sce- 
nario and our phantom based model, is the production 
of a huge amount of entropy in every cycle. Even in a 
"cyclic" scenario if one wants to address the usual cosmo- 
logical puzzles, such as flatness, homogeneity, etc, as well 
as produce cosmological perturbations with the correct 
amplitude required for galaxy formation, entropy pro- 
duction seems to be an inevitable requirement [261 ]. Like 
in inflation, in our scenario a huge amount of entropy 
is produced during reheating when most of the phantom 
energy is converted into radiation. This is essential to 
ensuring the cyclicity of energy density even though the 
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scale factor of the universe is monotonically increasing 4 . 

It is worth mentioning that in most cyclic models pro- 
posed the scale factor a(t) has a contracting-expanding 
behaviour, where both a bounce, near the Big Bang, and 
a turnaround, when the scale factor becomes large are 
needed. A brief and incomplete compilation of some of 
the papers that have proposed cyclic cosmologies can be 
found in 0, [l(| • Such cosmologies provide natural solu- 
tions to the flatness and horizon problems of standard 
Big Bang scenario. Some variants also avoid the issue 
of initial conditions, provided entropy produced during 
one cycle is not transferred to the next. In this case 
the cycles will not grow, i.e. become larger (Tolman's 
argument fill]) from one to the next, so we can no longer 
define a beginning of the universe. For instance the 
authors of [12| developed an ekpyrotic inspired cyclical 
model as an alternative to inflation, where a phase of 
slow contraction before the bang is responsible for gener- 
ating a nearly scale invariant spectrum of perturbations 
that seed the large scale structure formation. In con- 
trast, the standard inflationary scenario assumes a short 
phase that occurs after the big bang when the universe 
is rapidly expanding and nearly scale invariant perturba- 
tions are generated. 

The paper is organized as follows. In section |TT] we 
present our model of phantom fluid and discuss the cos- 
mology relevant for dark energy. In section Hnj we discuss 
how including interactions can lead to an inflationary 
space-time along with partially reheating the universe. 
In section HVl we provide a specific example where transi- 
tions from the phantom-inflation phase to radiation and 
vice-versa can be orchestrated giving us a cyclic model 
of the universe. In section [V] we discuss the different ob- 
servational constraints coming from inflation, Big Bang 
Nucleosynthesis (BBN), and dark energy experiments. 
Finally, we conclude with a summary of the scenario pre- 
sented and issues that needs to be addressed further. 



II. PHANTOM DARK ENERGY 

The purpose of this section is to implement a model 
for the phantom component that would drive the super- 
accelarated expansion. For now we will not include reg- 
ular matter nor radiation, but rather focus solely on the 
components necessary to obtain a phantom phase. One 
can possibly implement the cosmology sketched before 
in many different ways. Here we are going to realize the 
above picture using a ghost-like scalar field (with nega- 
tive kinetic energy) <j) coupled to a hidden matter sector 
which we denote by index "h" . There are two equivalent 
approaches to describe this type of coupling. One start- 



4 It may be possible to embed the scenario in phantom cyclic mod- 
els 0, |27H with actual phases of contraction, but we are not going 
to explore this possibility here. 



ing from an action, where we allow for a direct coupling 
term between the hidden sector and the scalar field. The 
second approach is to consider two fluids that can ex- 
change energy while maintaining the conservation of the 
total stress-energy tensor, as required by diffeomorphism 
invariance, although the individual stress-energy tensors 
are not conserved. 

Our system will be described by the following action: 



S 



M, 



R , (V0) 2 



(5) 



where Ch does not depend on <p 1 being the action de- 
scribing a perfect barotropic fluid. Here we work with 
a spatially flat FRW metric with signature (—,+ + +). 
Notice that the kinetic term for the scalar field comes 
with the 'wrong' sign, as appropriate for a ghost. In the 
phantom dominated phase, described here, the Hubble 
equation derived from the action above looks like 



H 2 = 



1 



1 



3A/ p 2 



Ph 



3A/ p 2 



(6) 



where we have assumed the field <p to be homogenous. 
Here a dot represents the derivative with respect to cos- 
mic time, t and ph denotes the bare energy density of the 
hidden sector, which is <f> independent. We will assume 
that it behaves like a perfect barotropic fluid, satisfying 
the continuity equation: 



p h + 3H(p h +p h )=0 



with an equation of state 



p h = uph 



(7) 



(8) 



One can also include a potential for the scalar field, and 
its effects are discussed briefly in the appendix, but for 
the purpose of illustration we are going to set it to zero. 

The interaction that we are going to consider between 
the hidden matter sector and <p is going to be very similar 
to the interactions considered in coupled quintessence (or 
interacting DE-DM) models [H US El] . From the action 
in equation ([5]) we get two additional equations of motion. 



Ph 



4> + 3Hj> = 2p h p(4>)/M v 
3H(l+Lu) Ph = 2p h 4>p((f>)/M p 



where we have defined 



(9) 
(10) 



(11) 



and we are going to assume hence forth that p(4>) always 
remains positive. As one can see, both the Klein-Gordon 
equation for cf> and the continuity equation for the hidden 
matter sector ph are augmented by interaction terms in 
the right hand side of the equations (|9ll0p . Although p in 
general depends on 0, to understand the phantom phase 
let us consider a constant p to begin with. It is easy to 
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check that the above interaction is consistent with the 
conservation of the total energy momentum tensor: 



It is as if the phantom scalar field is evolving under the 
influence of an "effective potential" given by 



Ptot + 3H(p tot + p tot ) = 



(12) 



where p tot = ph + p^. To see this we remind the readers 
that the Klein- Gordon equation can be recast as 



-2p h pcb 



3H( 



b) = -2p^{4>) 



(13) 



= PhW = PM^t-W' ( ^ ) 3(1+W) (18) 



An important thing to note is the +ve sign appearing in 
front of Ve' ff in equation (fTTj) because is a ghost field 
with negative kinetic energy. It is clear now that because 
of the peculiar properties of the phantom field, (f> actually 
rolls up the effective potential e 2 ^/ Mp . 

We have the following late time attractor power-law 
solutions: 



since the energy density and pressure for a phantom 
scalar field are given by 



P<j> = K 



(14) 



a(t) = a 



g iVX v — c p I 



Thus the source terms in the individual conservation 
equations (fTT))) and (fTB")) cancel each other. 

One can solve the hidden matter continuity equation 
exactly to find 



ph = PhoC(4>) — 



-3(l+w) 



(15) 



where we have chosen the convention that at a = a , 
C{4>) = 1 and ph = Pho- The a _3 ( 1+ ") dependence re- 
flects the usual dilution of the energy density of an ideal 
fluid with expansion. Depending upon the scale of en- 
ergy density relative to the mass of the hidden matter 
particles, they can cither behave as non-relativistic mat- 
ter (u) = 0) or like a relativistic species (to = 1/3) 5 . We 
will for most part consider a light degree of freedom, so 
that approximately it behaves like radiation. 

For the special case when p is a constant the Coupling 
function is given by 



C{4>) 



_ e 2 M (0-0o)/M p 



(16) 



with 



1-u 



4p 2 - 3/2(1 - lo 2 ) 



and p 



4/i 



4/z 2 - 3/2(1 - uj 2 ) 
(20) 

We have verified (see appendix [A| that these late time 
attractors are indeed stable 6 . 

In this phase the scalar field and the hidden matter 
are tightly coupled and evolve as a single fluid with an 
effective equation of state parameter 



Ph 



8 p 2 



p<j> + Ph 



3(1 



(21) 



The asymptotic value is attained during the late time at- 
tractor phase. In the phantom phase, expansion of the 

2 

scale factor is controlled by lu p , since a(t) ~ t 3(1+ "p' . 
These are analogues to the coupled quintessence solu- 
tions pjj [50, H3|. A detailed derivation of all results 
presented in this section can be found in the appendix 
[Al The crucial thing to note is that as long as 



Now, coming back to the evolution equations, we only 
need to solve the Hubble equation ((6]) and the Klein- 
Gordon equation, the latter simplifying to 



ZH<t> 



2ppho 



,2p(<t>-<t> )/M p 



-3(l+w) 



(17) 



>g(l- 2 ) 



(22) 



we have a phantom phase, i.e. 
for u> = 1/3 the last condition gives a constraint on p, 



< — 1. In particular 



(23) 



5 In this context we note that we have a choice in how we interpret 
the augmentation of the energy density with growth of fj>. One 
can either think of this growth as simply the increase in mass of 
the hidden matter particles if the mass depends on <f>, or creation 
of the hidden matter particles through its interactions with <f>, or 
a combination of the two. To keep things simple we are going to 
assume that the mass of the hidden matter particles remains a 
constant, but its number density increases. 



In passing we also note that in this phase most of the 
energy density is actually stored in the hidden sector; 



6 We have not investigated whether these solutions suffer from 
any hydrodynamic instabilities of the nature found in some in- 
teracting quintessence models [l9| , and we leave this for a future 
exercise. 
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one can check that the tracking ratio between scalar field 
and the hidden matter density is given by 



Ph 



K 
Ph 



8n 2 



3(1 - lu) 2 + 8p 2 



< 1 



(24) 



One can also calculate the amount by which cj> evolves 
during this phase. A straight forward calculation gives 



3(1 + w) 



hi 



Am 



(25) 



Here A m i n and A max represent the energy scale at which 
the phantom phase begins and ends respectively. 

In summary, in this section we have found under what 
conditions a phantom phase could be described by a sta- 
ble late time attractor for ghost-like fields coupled expo- 
nentially to a perfect fluid. The purely phantom sector 
we have studied in this section is problematic as an in- 
flationary model. For example, for a constant (purely) 
phantom equation of state u> p < —1, typically one would 
obtain a blue spectrum (see equation (|59|) ) which would 
be inconsistent with the WMAP data. This conclusion 
is of course valid with the assumption that primordial 
perturbations are generated mostly during the phantom 
phase. Thus unless p is fine tuned to be very close to the 
critical value (which gives rise to u) p = — 1), we will not 
be able to reproduce the inflationary near scale-invariant 
spectrum. Thus far we have considered the phantom sec- 
tor alone; in the next section we include interactions with 
the standard model which ameliorate some of the prob- 
lems of a phantom sector alone. 



III. "PARTIAL REHEATING" AND LATE TIME 
DESITTER PHASE 

In the previous section we realized the phantom phase 
through an interacting phantom scalar field and hidden 
matter sector. In the absence of any new physics this 
phase is going to last till the Big Rip singularity, as is well 
known in phantom cosmology. In order for the next cycle 
to begin we need to first find a "reheating" mechanism 
which converts most of the phantom energy density to 
radiation. 

What we find, quite remarkably, is that once we in- 
clude interaction between the hidden matter sector and 
Standard Model particles (namely we allow the hidden 
sector particles to be converted to light degrees of free- 
dom of the standard model), which generically exist, it 
naturally ameliorates the phantom like acceleration to a 
near exponential inflationary expansion. 

As we will see, such interactions begin the process of 
reheating the universe by producing a radiation bath. 
Unfortunately the interactions don't provide us with a 
graceful exit from the inflationary phase, but we will dis- 
cuss how this issue can also be addressed in the next 
section. 



To understand how interactions effect the cosmological 
evolution we will use Boltzmann equations in the follow- 
ing form: 

p h + 3H(l + w)p h = -p h r + 2p h j>-£- (26) 

Mp 



4> + ZH<i> = 2 Ph 
p-y + AHp 1 = p h T 
along with the Hubble equation 



M„ 



H- 



1 



3M p 2 



Ph 



Pi 



(27) 
(28) 

(29) 



Here by p 1 we include all the light degrees of freedom 
which do not couple to the phantom field, and T is the 
annihilation rate of the hidden matter particles into all 
these other light degrees of freedom. We have ignored the 
inverse process of creation of the hidden matter particles 
from the rest of the matter under the assumption that the 
equilibrium density of the hidden matter sector is small 
compared to normal radiation. Now, the annihilation 
rate, per hidden sector particle, is given in general by 



hh — >77 



rih < <j\v\ > 



hh — >77 



(30) 



where < a\v\ > is the average over all initial and final 
states of the differential cross-section times the relative 
velocities of the annihilating particles. Usually one is 
used to consider the opposite process while trying to de- 
termine when a given species freezes out. In the latter 
case since the photons are in thermal equilibrium, one 
can use thermal distribution functions to compute the 
"thermally averaged" < a\v \ >. However, in our case in 
order to compute < a\v\ > we would need information re- 
garding the velocity distribution of the produced hidden 
particles from (f>. In the absence of any micro-physical 
theory of such an interaction, for the purpose of illustra- 
tion, here we are simply going to assume that < <r|i;| > 
is a constant set by the details of the interaction, so that 
the interaction rate per hidden sector particle is given by 



Ph 



(31) 



Here m is an energy scale we introduce as a free param- 
eter. More generally one expects T to go as some power 
law, r ~ p^, the power being determined by the micro- 
physics, but most of our results and conclusions should 
hold qualitatively as long as A > 1. 

It is easy to check that the above set of equations (|264 - 
[28]) have an asymptotic de Sitter late time attractor so- 
lution where all the energy densities and the Hubble pa- 
rameter tend to a constant. Defining the following di- 
mcnsionless variables, 



Ph 



M 2 H 2 



and 



r 

H 



Ph 
m 3 H 
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we have: 



27w - 9 + 3^9(3^ - l) 2 + 192fi 2 



3w - 9 1 



Pi 



2 

ZphP 
3H 

Pi 
AHm 3 



-v/9(3c-l)2 

1\jbrr?r\ 
3M p 
m 7] 



192^ 2 



4M 2 £ 



(32) 
(33) 
(34) 

(35) 



As a consistency check we have evolved these equations 
numerically, sec Fig. 1 . and have verified that the asymp- 
totic values are exactly the ones predicted by the above 
set of equations. Let us make a few observations. Since 
we are specializing to w = 1/3, let us look at the asymp- 
totic value of the energy densities in this case: 



Pi 



Ph 



Pd, 



m 

9M 2 

16 V3m m 

9M| 

32 fj, 2 to 6 
9Mf 



V3/i- 1 



(36) 



Note that the solutions above are consistent only if 
fi > -^j, but we know this is also a requirement for the 

existence of a phantom phase 7 which will eventually set- 
tic to this dcSittcr attractor. In particular, wc find that 
for ^> 1, which is the case we will eventually be focusing 
on, radiation dominates over the hidden sector: 



P 1 ~ 4V3> 



< 1 for n > 1 



(37) 



In other words, as \i increases, the conversion from hid- 
den matter to radiation becomes more efficient, and in 
particular radiation can easily dominate the total energy 
density. However, this does not mean that we can en- 
ter into the usual radiation dominated epoch because 
the equation of state for all the energy densities essen- 
tially approach tu p — — 1, as all the energy densities ap- 
proach the constant asymptotic values in equation (|36p . 
Physically, the energy density that the phantom field was 
pumping into the hidden matter sector is now transferred 
to radiation through its interaction with the hidden sec- 
tor, in such a way that we approach a deSitter universe. 
This phase however is good for inflationary cosmology 
as fluctuations produced during this de Sitter phase are 
expected to give rise to a near scale-invariant spectrum 
(allowed, although not favored, by the WMAP data). In 
order for the amplitude of the fluctuations to be consis- 
tent with observations we require 



p 7 



10~ 12 M?, 



1(T 2 M B 



(38) 



We will return to the constraint coming from the num- 
ber of efoldings required for a successful inflationary 
paradigm later. 

In this section we have found that the big rip singu- 
larity could be avoided in this model if the hidden sector 
particles are converted to light degrees of freedom of the 
Standard Model. As the energy densities approach their 
asymptotic values the universe will enter in a deSitter 
phase. Of the three components, radiation is dominant, 
yet the universe is inflating. This is due to the interplay 
between the coupling of the hidden sector to the phan- 
tom field and light degrees of freedom of the Standard 
Model which leads to a state where all the energy densi- 
ties approach a constant value. 



IV. TRANSITION TO RADIATION AND 
CYCLICITY 



In the above section we saw that interactions between 
hidden matter and radiation can ameliorate the phantom 
like acceleration to exponential inflation but cannot pro- 
vide a graceful exit from the dcSittcr inflationary phase. 
So far we have assumed that fi remains a constant lead- 
ing to an exponential coupling between the hidden sector 
and the phantom field as in equation Q16p . If fx(4>) is not 
a constant a graceful exit becomes possible. Here we ex- 
plore the case when /J.(<f>) is periodic. For simplicity we 
will just assume a step function for /x(0), where 



/i = jU p 3> 1 for < <f> < 4)r 
fj, = fi r <C 1 for 4>r < <f> < 0o 



(39) 
(40) 



sec equation (1 23 I t 



and then the pattern repeats itself. The first phase, when 
/j, ^> 1, reproduces the phantom phase discussed in the 
previous section, leading eventually to the late time de 
Sitter like attractor evolution. However, now this infla- 
tionary phase ends once cf> reaches the transition value 
4>r. We are also going to assume that u = 1/3 in this 
scenario. We remind the readers that ui describes the 
hidden sector, as in equation (JH). As we will see shortly, 
the periodicity in (i will ensure that we enter the stan- 
dard radiation dominated era which lasts till rolls to 
</>o- 

The evolution of the universe during one cycle can be 
described in our model using three phases, as indicated in 
Fig. 1. Phase I provides the "Reheating" from inflation- 
ary expansion. Phase II describes a standard radiation 
dominated era. Then a short phantom phase III A ensues, 
during which the energy is driven from the meV to the 
GUT scale. This is followed by a deSitter inflationary 
phase IIIB. After one such cycle is complete, the next 
one begins, again cycling through Phases I, II, IIIA, and 
IIIB successively. 
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A. Phase I: Reheating 

The reheating process is most easily understood when 
p r = (defined in equation (j40|) ). so let us first focus on 
this simple case. As we discussed in the previous section, 
as long as p is large, although radiation is the dominant 
energy density, the deSitter phase continues. However, 
if and when p sharply falls to zero, the phantom phase 
indeed ends. Two things happen. Firstly, since initially 
r is comparable to H, as it can be seen from equation 
(|33p . there is rapid conversion of the hidden matter to 
radiation, but the hidden matter sector now no longer 
gets replenished by the scalar field. Secondly, the driving 
term in the right hand side of the Klein-Gordon equa- 
tion for the phantom scalar field (J9j) is now absent and 
as a result (j> slows due to Hubble damping, p^ ~ a~ 6 
and eventually comes to a halt. At the beginning of the 
reheating phase the energy densities of radiation, hidden 
matter and <p are approximately given by the asymptotic 
values of the late time de Sitter phase (f3"o| . As we had 
pointed out before, for /i> 1 radiation is the dominant 
energy density component in this asymptotic phase. The 
"reheating" phase further ensures that radiation contin- 
ues to dominate the energy density. If one tracks the ratio 
of the energy densities between hidden matter and radia- 
tion, J?, then it starts with the asymptotic value given by 
equation (|37[) . decreases rapidly during conversion, and 
then approaches a constant, ^ m ± n , once the conversion 
ends. Note, since slows down, ph redshifts almost as 
radiation and therefore maintains an approximately con- 
stant tracking ratio approaching & m ± n . In appendix [Bl we 
calculated this asymptotic ratio to be 

1 /3 

^mm ~ 7, — for Mr = , p p > 1 (41) 
o p p 

Qualitatively, it turns out that one can distinguish two 
different regimes in the reheating phase depicted in Fig.l 
as phases LA and LB. Numerically we found that even 
after p — ► 0, it takes a while for the radiation energy 
density to start decreasing substantially. The reason is 
somewhat technical and the reader is referred to the ap- 
pendix [B] for details. Intuitively, the main reason is that 
initially the Hubble damping of radiation is compensated 
by the hidden matter decays into radiation, AHp^ ~ p/iL. 
Since radiation is the dominant component of the energy 
density, this in turn leads to H being approximately con- 
stant, as can be seen in phase LA of Fig.l. Once ph 
decreases appreciably so that phX <C 4Hp 7 , the radia- 
tion energy density starts to decrease appreciably and 
therefore so docs the Hubble rate. This is depicted in 
phase LB of Fig.l. At some point, the conversion from 
hidden matter to radiation effectively stops, marking the 
end of the reheating phase. 

Although the above discussion has been for p r = 0, 
we note that for a non-zero but small p r , the reheating 
phases (LA, B) follow basically the same pattern as in the 



p r = case. The asymptotic tracking ratio 8 between 
hidden matter and radiation receives a slight correction: 

^ 1 V3(3 + 8u r ) „ , 

.^min ~ ^ 2 ' for p r < 1 , p p > 1 (42) 

to Pp 

We have checked this numerically, and some of the more 
technical details are discussed in Appendix [B] 

B. Phase II: "Standard" Radiation Domination 

After the reheating phase, since hidden matter is no 
longer converted into radiation, the latter starts to evolve 
as a(t)~ 4 , and consequently H ~ l/2i as in the standard 
radiation dominated era. In the meantime <j) continues 
to be Hubble damped. Once the scalar field effectively 
stops evolving, the hidden matter starts redshifting as 
radiation and thus settles down to its constant tracking 
ratio given in equation (|4ip . In particular, one can see 
from (|4"2"|) that for sufficiently large values of p p , this ratio 
can be quite small and easily satisfy constraints coming 
from BBN and CMB. BBN/CMB only constrains the 
abundance of dark radiation component to be less than 
around 10% [Hj]. In Fig.l we refer to the phase when <fi is 
being Hubble damped as phase II. A, and the subsequent 
radiative phase as phase II. B. 

If p r is precisely zero, then the radiation dominated 
phase IIA can continue forever because <p will effectively 
come to a halt, and unless the value of 4>o is fine-tuned, (f> 
will never make it to the next large p phase. As a result 
the next phantom phase will not begin and the cyclic 
picture cannot be sustained. This is why a small but 
non-zero value of p r is essential to maintaining cyclicity 
without having to resort to unnatural fine-tuning. 

For a non-zero but small p r , the reheating phases 
(LA, B) and the Hubble damping phase (II. A) follow very 
much the same pattern as in the p r = case, as discussed 
at the end of the previous section IIV Al The main dif- 
ference when p r 0, as compared to the p r = case, 
appears in the radiative phase II. B. This is because the 
non-zero driving term in the Klein-Gordon equation now 
ensures that instead of coming to a halt, <p now tracks ra- 
diation. After the initial phase of Hubble damping, the 
driving term on the right hand side of the KG equation 
catches up with the Hubble damping term. At this point 
the scalar field enters a phase where its energy density 
approximately tracks that of radiation. This can be seen 
in Fig.l where we have divided the radiation phase into 
two parts. IIA refers to the regime when the scalar en- 
ergy density is still being Hubble damped, while phase 



8 Unlike in the fi r = case where the tracking ratio keeps decreas- 
ing and asymptotically approaches ,^ m i n , when /i r ^ 0, .S^min is 
actually a minimum of the tracking ratio that is attained. Since 
never comes to a halt, the ratio does increase from its minimum 
value of ,^ m in, but this increase is very slight. 



II. B refers to the tracking phase where both the hidden 
matter and the scalar energy densities are tracking that 
of radiation. 

This tracking behavior can be approximately obtained 
as follows: From the KG equation we have 

3H(j> = w — w 6M mia M p p r H (43) 

IVlp Mp 

In the above we have ignored the e 2Mr * dependence 
of ph, as <f> is rolling very slowly, and p, r <C 1. We will 
perform a consistency check later. Also, we have ignored 
the contributions to the energy density coming from the 
hidden matter and the scalar field as compared to normal 
radiation. Again, this is justified as 7?. m i n <C 1 and (f> is 
rolling slowly. Choosing the ansatz 

<l> = M p p r ln(£)^j>=^f L (44) 
we find that (|13")) can indeed be satisfied provided 

Pr = .^minMr (45) 

In the above analysis we have used the fact that in a 
radiation dominated universe H ~ l/2t. In particular 
our analysis tells us that the tracking ratio between the 
kinetic energy of <fr and radiation is indeed very small 

^ = -liPr^nf (46) 

justifying our earlier assumption. In the next subsec- 
tion we will also see that during this phase <p evolves 
rather slowly, so that C(4>) changes only by an 0(1) fac- 
tor ensuring that the hidden matter indeed behaves as 
radiation to a very good approximation. 

This radiation dominated tracking era lasts till <j) 
reaches <po an d rolls over to the large p region. The next 
phase of phantom domination, phase III. A, can now be- 
gin. 

C. Cyclicity 

To better understand the transition from one cycle to 
the next let us discuss the various phases we observe in 
Fig. 1, where we plot a numerical solution for the en- 
ergy densities of the three components from one deSittcr 
inflationary phase to the next. The plot does not corre- 
spond to realistic values for A max or A m i n , but captures all 
the essential features of the different phases. The plot 
starts (extreme left) at t = and <fi = <pn, corresponding 
to the beginning of the reheating phase LA. We have 
estimated in Appendix [B] equation (|B8|) how long (£ia) 
the LA sub-phase lasts 

/ 3M 4 A ' 
i 777O77J 1 
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FIG. 1: Numerical solutions for the energy densities as we 
complete a cycle from a deStiter phase back to it. Here we 
have chosen: u) = |, p v = 5, p r = .1, m = 10~ 2 M P and we 
have set M p to one. Note the six distinct phases: LA and LB 
corresponding to reheating; II. A and II. B corresponding to a 
radiation dominated universe; III.A and III.B corresponding 
to the phantom and dS phase respectively. In order to make 
all phases clearly distinct we chose the transition set the min- 
imum energy density at around 10 -45 instead of the realistic 
meV 4 

During this phase the field <fi evolves approximately a 
distance of: 

Aj> 1A = <p dS die" 3 "-* = 4g-[l- e- 3H ^] 
Jo oH ds 

(48) 

Above we have used: 

HhA) = <Pdse- 3HdstlA (49) 

where 4>dS represents the asymptotic value in the deSittcr 
phase, given by ([34]) and Hds is the Hubble rate during 
the inflationary phase and can be solved for from the 
definitions of £ and 77. 

Phase LA is followed by phase LB, where although the 
conversion from hidden matter to radiation takes place, 
the Hubble rate starts to decrease appreciably as well. 
As soon as the conversion is no longer efficient we enter 
a regime, phase II. A, where radiation starts to rcdshift 
as a~ 4 marking the beginning of a Standard radiation 
dominated era. In Phase II. A radiation and the hidden 
sector energies approximately track each other while the 
field 4> is still being Hubble damped. This phase ends 
when the scalar field is no longer Hubble damped and 
starts to track radiation as well, a phase we refer to as 
II.B. Next we estimate the time <2A, when this transition 
from the Hubble damping phase II. A to tracking phase 
II.B occurs. It can be defined as the time when the Hub- 
ble damping term in the l.h.s. of equation ([27)1 is equal 
to the coupling term on the r.h.s. Under the assumption 
that we are in a radiation dominated phase and that the 
hidden sector energy density tracks the radiation energy 
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density with the ratio ^min we get 



2*2,4 



M p 

t\ A 4> 2 {t 1 A ) 



2 

min 



(50) 



This last equation can be rewritten in terms of our pa- 
rameters in the following form: 



t_XA 



1 p 2 r (3 + 8p r ) 2 r, 



fiH d sti 



144 



(51) 



With u> = 1/3 and in the limit p, p 3> 1 the exponent 
&Hdst\A becomes v^, so we see, as expected, that t^A ~> 

tlA- 

We can compute the distance the field 4> evolves during 
the phases LB and II. A since the Hubble friction term is 
dominant during this time. 



dtcp = 2<p(t 1A )t 




1 




(52) 



In the next phase, II. B, as one can see from the plot, 
all the energy components are tracking each other. Since 
this phase will lasts until almost 'today' and it began at 
t2A we have: 



A0 : 



2B 



M„ 



20L 



(p 7 (W) 4 
meV 



120^ ln /i r (53) 



where A<f> 2 B is the distance the scalar field evolves during 
phase II. B. In particular we note that this means 



(54) 



Since & m ±n is a small number, for sufficiently small val- 
ues of fi r it is easy to see that the exponential will only 
contribute to an 0(1) factor to the energy density of the 
hidden matter sector. In other words unlike the fi r = 
case, although the ratio S% does not monotonically de- 
crease during the radiation phase, but rather starts to 
increase as <fi evolves, this increase is very slow. This 
justifies our earlier assumption of hidden matter approx- 
imately behaving as radiation. 

The advantage of having a non-zero p r is that it keeps 
the scalar field rolling, albeit slowly, ensuring passage to 
the next phantom phase when </> reaches (fig. Therefore, 
no fine-tuning is involved in restarting the phantom era. 
To see this observe that 



2B 



bo - 4>r) - (A<fn A + A0 



IB+2A, 



(55) 



and crucially for realistic values of the parameters the 
three different A</>'s are of the same order of magni- 
tude. If for instance, it turned out that A<f>2B "C 



A<f>iA + A(f>iB+2A, that would have meant fine tuning 
the range O — <t>R to cancel A$ya + A<f)i B+2 A to very 
high precision. In fact, this is what one has to do as 
fi r — ► 0. 

A related nice feature of the model is that the expo- 
nential hierarchy between the scales of inflation and dark 
energy is rather easy to arrange. To see this let us try to 
obtain A min in terms of 4>q — <j)R. By re-arranging equation 
(|53| one finds that the energy density at the beginning 
of the phantom phase is given by 9 



A m 



1CT 3 M P exp 



A m 



: exp 



A^, 



2B 



2&minPr M p 



b R - A<j) 1A - A<j) 1B +2A 



(56) 



Since & m in,Pr are small numbers as compared to 
A0/M p 's, it is easy to arrange the exponential suppres- 
sion of po as compared to the GUT scale. Conversely, we 
can reformulate the "smallncss" problem associated with 
dark energy in terms of four parameters: (i p , jx r , 4>r, (po, 
which as we shall see shortly have values that range from 
0(1) to O(10 2 ). In the next section (below equation ([62])) 
we will provide specific numerical examples which make 
this more evident. 

Finally, we come to the phantom phases III.A,B. The 
transition to the phantom phase, which is supposed to 
be happening during the present cosmological epoch, oc- 
curs once (j> reaches <po and [i transitions suddenly to its 
high value, fi p . The hidden matter + scalar field en- 
ergy densities catch up with radiation, start evolving as 
a phantom fluid with equation of state given by (f2"Tj) , and 
come to dominate the universe. The radiation keeps get- 
ting diluted as a~ 4 during this phase which we refer to 
as phase III. A. Since the increase in energy density in 
the phantom phase occurs at a very short time scale, all 
the features of this phase are not very discernable in the 
log-log plot in Fig.l, but we have checked them numeri- 
cally As we have seen before, the annihilation term will 
ultimately cause this behavior to transition to a de Sitter 
phase, when r ~ H or ph ~ m 6 /AI p . At this point we 
enter, what we call phase III.B, when all the energy densi- 
ties become comparable and then tend towards constant 
values, leading to an asymptotic de Sitter space-time. It 
is during this phase we expect to generate a scale invari- 
ant spectrum of perturbations. This phase will last until 
the field <p has evolved a distance of tfrg + 4>r, when we 
re-start the cycle. 



9 Above we have used a value for p^ 4 (t2A) of 10 — 3 Mp, slightly 
overestimating it. The actual value is typically lower than this, 
because of the additional Hubble dilution of the energy density 
of radiation from t%A to t?.A- For instance, with the parameters 
used in the numerical solution for Fie 



10- 5 M B 



1 we have p^ 4 (t2A) 
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V. NUMBERS AND CONSTRAINTS 

So far we have provided general constraints coming 
from different observations on the couplings and the 
scales. For the purpose of illustration let us provide some 
typical values which conform to these constraints and in 
the process we will also be able to understand the dif- 
ferent phases of evolution better. Let us start with the 
inflationary phase. As already discussed, to obtain the 
correct amplitude of inflationary fluctuations we will take 
the reheating to occur at approximately GUT scale en- 
ergy densities which implies 



10" 2 M„ 



(57) 



During the course of 60 c-foldings 10 , we find, using 
equations p3|34|) . that the field <f) evolves a distance of: 



A<jy 3B = 40A V P £ 



(58) 



As a prototype example, for /i p = 4 this gives us Afos ~ 
200. This only gives us a constraint on the range when /i 
is large. For instance, for n p = A the number of cfoldings 
during the phantom phase turns out to be (using equa- 
tion ([23]) ) A<p3A ~ 30. To be consistent with inflation 
we must have <j)R > 200 + 30 = 230. Given the fact that 
these field values are transPlanckian, it is possible that 
higher order terms in the Lagrangian should not be ne- 
glected; we proceed here with the assumption that our 
starting point is sensible nonetheless. 

Depending on the number of c-foldings during the de 
Sitter phase we find two distinct cases for the spectrum 
of fluctuations. Case I. If the range of <fi is such that 
one gets more than 60 e-foldings of de Sitter, then the 
CMB fluctuations are scale invariant. As we noted in 
the introduction, n s f=a 1 is still consistent with obser- 
vations 0] once one allows the possibility of running of 
the tilt and/or tensor modes. If the range in <fi is such 
that we have only around 60 e-foldings then the CMB 
fluctuations at large scales can show a transition from a 
blue (when the phantom phase will be operating and the 
reheating mechanism hasn't kicked in completely) to a 
scale invariant spectrum. This could be an interesting 
and rather unique signature of the model. Case II. If the 
number of cfoldings in the dS phase is shorter than 60 e- 
foldings, the fluctuations that we are observing in the sky 
must have been generated in the phantom phase. This 
gives a rather stringent constraint on how far below —1 
the phantom equation of state uj p can be. For a "super- 
inflationary" space-time sourced by phantom fluid, the 
spectral tilt is expected to be given by 



'Is 



6(1 + cj p ) 
1 + 3w„ 



> 



(59) 



10 We choose 60 as a generic typical example for GUT scale infla- 
tion; the actual number of e-foldings might be lower or larger 
than this. 



implying a blue spectrum. Therefore to be consistent 
with the observations uj p has to be very close to — 1 . For 
instance, if we include tensor modes, according to at 
the 2cr level we find 



Vs < 1.01 



> -1.01 



(60) 



implying fi p < .581. This bound is very restrictive, re- 
membering that in order to have a phantom phase we 
need \i v > 1/V3 « 0.577. In the rest of the section we 
are not going to discuss this possibility any further and 
concentrate on Case I with fj, p 3> 1 which seems more 
attractive. 

Let's next look at the constraint coming from BBN 
and WMAP on the amount of dark radiation [T^ |. To 
be consistent with the data the amount of dark radiation 
has to be limited to within 10% of ordinary radiation. 
This essentially imposes a constraint on /i p , but a rather 
weak one 



< 0.1 



/'/< 



> 0.17 



(61) 



which is easily satisfied. 

What about the range of A<j) during radiation domina- 
tion? This of course depends on the value of [i r . Just 
to get an idea, we find for /i p = 4, \i r = 0.2, we need 

A(f) = A(f)i A + A(/)ib+2A + A(/) 2 B = (f>0 — 4>R ~ 

- 0.70 + 2.16 + 0.66 (62) 

As one can see, all the values in the above equation are of 
the same order of magnitude and therefore no fine tuning 
seems to be involved. 

Thus the picture that emerges from the above esti- 
mates is that for the scenario to work we need a rela- 
tively longer phases in <f> when /i is large (|58|) . followed 
by shorter phases when [i is small. However, even the dis- 
crepancy in Acj) between equations (|58[) and (|62|) is only 
a few orders of magnitude. Similarly, fi r and fi p (in the 
above example) differ again by only a few orders of mag- 
nitude. It is these numbers (the values of A^ and the 
two values of fj,) that determine the heirarchy between 
today's meV scale and the GUT scale of inflation via 
equation |56|) . So we see that, indeed, little fine-tuning 
is required to explain these disparate mass scales. 

Finally, let us come to constraints from dark energy. 
The main constraint comes from the equation of state 
parameter. A combined 2er bound from CMB+BAO+SN 
is given by —0.88 > uj p > —1.14 0. This provides a 
constraint on \i p 



H P < 0.61 



(63) 



This may seem too small, but we realize that currently 
we are undergoing the phase transition from small /i to 
a large fi region, and therefore it is easy to arrange that 
the "current" value of \i satisfies inequality (|6lS|) and has 
not yet reached the constant maximum value n p . Equiv- 
alently, the equation of state parameter today has not 
yet reached its late time phantom phase value given by 
equation (|2"Tj) . 
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VI. CONCLUSIONS 

In this paper we have studied a cyclical model of the 
universe where the energy density cycles between a min- 
imum value, typically of the order of meV 4 and a maxi- 
mum value roughly set by the GUT scale. This effectively 
provides a connection between the current accelerated ex- 
pansion we observe today and the inflationary era in the 
past. The scale factor continues to grow from one cy- 
cle to the next (there is no "turnaround"). In order to 
achieve this model we postulated the existence of some 
'hidden' sector matter coupled to a ghost like scalar field. 
This mechanism is responsible for a super-accclcratcd 
phantom expansion. Allowing for hidden sector parti- 
cles to be converted to light degrees of freedom of the 
standard model ameliorates the phantom behavior, effec- 
tively transitioning to a dcSittcr like expansion, therefore 
avoiding the Big- Rip singularity. Although dominated by 
radiation, in this phase all the energy densities remain 
constant. Therefore, if most of the cosmological pertur- 
bations are generated during this exponential inflation- 
ary era the spectrum is expected to be scale invariant. 
Even if not favored by the data, this is still a possibil- 
ity. Another possibility would be to have some of the 
fluctuations generated during the phantom phase which 
will show up as a blue tilt in the spectrum. In that case 
we will see a running of the tilt which could be a unique 
possible signature of the model. 

In order to achieve the cyclic behavior we have postu- 
lated that the coupling between the ghost scalar field and 
the hidden sector is constant piecewise. This procedure 
might seem ad-hoc, but it is just the simplest possibil- 
ity. We found that no fine tuning seems to be involved 
when requiring to have a long enough radiation/matter 
dominated phase. It is also worth mentioning that the 
"smallnesss" problem associated with dark energy is cir- 
cumvented. The only parameter we need in order to de- 
scribe the current acceleration is the coupling between 
the hidden sector and the ghost field, and it has a value 
not much greater than one. 
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APPENDIX A: STABILITY OF SCALING 
SOLUTIONS IN PHANTOM COUPLED MODELS 

We will look at ghost fields coupled to radiation or 
matter via the term p 1 = p 7 e 2/ ^ and study the stability 
of the critical points, generalising the analysis done in 
[23d . In this way we will find the stable attractors and 
conditions necessary to obtain them. Here the subscript 



7 refers to the adiabatic index of the fluid the phantom 
field couples to, defined as 7 = 1 + u>h- Throughout the 
main body of the paper we have used ph, but for nota- 
tional convenience here we will replace it by p 7 . It will 
be assumed that the phantom field has a self-interaction 
potential of the following forms: V(<f>) = Vbe -2 "^ and 
V{4>) — —Voe 2a ^. We have not used a potential in deriv- 
ing our main results in the paper, but we keep it here for 
generality and further reference. 
For 



V{4>) = V e- 
we have the following equations: 
4> + 'AH4> = VeffW 



H 

H 2 



1 



(Al) 



(A2) 



-y +F(0)+p 7 e 



where V e ff = V(cf>) + p 7 e 2/i ^. Although we are treating 
here only two fluids, this case is relevant during the phan- 
tom phase, when the energy density of any (third) com- 
ponent not coupled to the phantom field, such as regular 
matter, will quickly become sub-dominant. Introducing 
the variables 



and 



y 



V6H 



(A3) 



(A4) 



and using T) = log a(t) as the independent variable, in- 
stead of the cosmological time t, we can rewrite the equa- 
tions in the form of an autonomous system supplemented 
with the Hubble constraint. 



-3x + V6[n{l 



V 2 ) - ay 2 ] + 



2x z 



V = 



3 

+ -x [^(l + x 2 -y 
~V6axy+^y[j{l + x 2 -y 2 )-2x 2 ] (A5) 



2 2 
y - x 



In order to find the critical points one has to set the 
r.h.s of the first two equations to zero and solve for x and 
y. We find five distinct solutions, 



,v) 



(i 

H 

V 3 -2+7 
( 1 

\ 4 a+fi 



0) 
0) 
0) 



1 Vl^Sfia+Sfj, 2 - 
a+fj. 



■67-I-37 2 



r%/3V3 + 2a 2 ) 



I. 
II. 

III. 

IV. 
V. 
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Note that the first two are nonphysical so we will no 
longer consider them. Next we compute the fractional 
densities for the 'phantom' field and its adiabatic con- 
stant near the five critical points. The Hubble constraint 
will then enforce additional existence conditions. 



III. 



^0 = y 2 - x 2 



"3 (-2+7)* 
1 — 4/iQ — 4/i"+37 j-y 



l 



(a+M) 2 



(A6) 



V. 



One solution is completely dominated by the scalar field 
4> and the other two exhibit a scaling behaviour. For the 
first of those the Hubble constraint will not give any ad- 
ditional inequalities in the parameter space since Sl^ni) 
is clearly negative, hence less than 1. Here the roman 
numeral subscript refers to all the solutions, including 
the non-physical ones, i.e (/) corresponds to the solutions 
with [x, y) = (i, 0) and so on. Also CImiv) < 1 is trivially 
satisfied for positive parameters. So the only nontrivial 
existence constraint comes from imposing reality of the 
yiv solution: 



juOu + a) > -7(2-7) 



(A7) 



Let us now look at the 'effective' values of the adiabatic 
constant 7 = 1 + w for the ghost field near the critical 
points. For completeness we will keep even the nonphys- 
ical solutions. 



J<t> 



1 + i^ 

P<t> 



2x 2 



2 - V{4>) x 2 + y 2 



(A8) 



With this definition we get, 



2x 2 



70 



-y 2 + x 2 



' 2 
2 

= \ 2 



37' 

—iua— 4^+37 

■±a 2 
3 



/. 
II. 
III. 

IV. 
V. 



(A9) 



The total, or 'effective' DE equation of state parameter 
that drives the expansion can be defined since we have 



scaling solutions. 



Utot — 



V<p+P~t 



(A10) 



P<t> + P~i 
Some algebra leads to 

utot = fy(7*-l) + (l-fy)(7-l) ( Al1 ) 

For the three physical solutions it leads to the following 
values: 



1 3w h (ui h -l)+8n J 

3 



Wtot 



U3h — 1 

= <^ 1 3(l+c^) 2 +2(^+a)(^ h q- M ) 
2 (5+7T? 

w 



1 



III. 

IV. 

v. 



(A12) 



Next we will study the stability of the relevant critical 
points, ignoring the non-physical first two solutions. The 
technique is the following. One expands around the crit- 
ical solution setting x = x c + u and y = y c + v into 
(|A5[) and then keep only linear terms. In order to have 
a stable solution the eigenvalues of the matrix describing 
the linearized system must have negative real parts. For 
solution III. one finds the following equations 

, 1 - 37 2 + 12 7 - 12 



2 -2 + 7 

1 (-37 2 - 8/j.a - 8p? + 67) 

2 -2 + 7 1 



(A13) 



For both radiation (7 = 5) and dust (7 = 1) the coeffi- 
cient u is negative independent of the value of /i. From 
the second equation one gets that the node is stable when: 



37(2-7) < 8/i(a + /i) 



(A14) 



Otherwise we have a saddle point. It is interesting to 
notice that this is exactly the same condition we got in 
(|A7|) for the existence of the fourth solution. 

Let us move on to the stability of solution IV, for which 
the linearized system becomes: 



, 1 24a 2 - 24 7 ^ 2 + 24/r 2 + 48/xa + I87 2 - 36ja[i - 9j 3 - 12 7 a 2 

U = ; rr; U 

8 (p + a) 2 

1 V3y/8fia + Sf? - 67 + 37 2 (16/ia + 37 2 + 8fi 2 + 8a 2 ) 

2 v 



8 (a + n) 

1 V^j8na + 8fi 2 - 67 + 37 2 (-4a 2 + 37 2 - A^a - 67) 3 M8ua + 8^ 2 - 67 + 37 2 ) / 4 lrN 

; rr; U ; V (A15 

8 (a + /i) 2 8 (a + /i) 2 y ' 



The eigenvalues are: where B is the following combination of the parameters: 

() 4a + ^V ' p J W(36 + 8l7 2 - 180 7 ) + 367(7 -a/i + 4/i 2 ) 

e 2 f 4 ) = Uai- 6a- 6u- B?) (A16) 

4 a + /j, \ ) 
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The study of stability of the solutions is quite compli- 
cated here but there is a range of parameters for which 
the real part of the two eigenvalues is negative. Inde- 
pendent of the value of B one necessary condition for 
stability here is 



2(1 + ^) < 7 
a 



(A17) 



For positive parameters, and if 7 is cither 1 or | this in- 
equality cannot be satisfied; thus the fourth attractor is 
unstable if the hidden sector is comprised of some com- 
ponent that behaves like matter or radiation. 

Let us look now at the stability conditions for the crit- 
ical point labeled by V. The linearized autonomous sys- 
tem becomes in this case: 



u! = (-3 + 2a 2 (7-3)-4/ia)M + 

+v / 2v / 3 + 2a 2 (2/i - a 7 + 2a)v 
v' = V2\/3 + 2a 2 {aj - a)u - 7(3 - 2a 2 )i 



(A18) 



The eigenvalues read — (3 + 2a 2 ) and —(4a 2 + 4fia + 37) 
which are clearly both negative and real so one has a sta- 
ble node as a late time attractor. Also remember that 
this solution does not have a scaling behavior, since in 
this case the energy density is dominated by the phan- 
tom. In conclusion we found two stable late time attrac- 
tors, one of which exhibits a scaling behaviour. Notice 
also that this scaling solution corresponds to a critical 
point where y = 0, so in effect it is equivalent to the case 
where the potential is actually zero. 

Next let us look at the tracking ratio in the two cases of 
interest, namely the third and fourth critical points. The 
tracking ratio is defined to be the ratio of relic densities 
of the two components. 



Sid 



1 -si, 



(A19) 



With this definition we find, using (|A6 



3 (-2+7)' (i+s 

1 -4^q-4/j 2 +37 



(/x+a) 2 1 + 



III. 



IV. 



Since we have established that wc will be using the 
late time attractor described by critical point III. as our 
model for the phantom phase, let us actually see under 
what conditions we cross the 'phantom divide' and if this 
phase will be stable. Comparing u> tot given by (|Al2j) the 
following restriction on fi: 



2 1 | for matter hidden sector 
' I for radiation hidden sector 



(A20) 



Looking back at (|A14[) and taking a = as appropri- 
ate for the case of no potential one can check that the 
stability conditions are actually identical with the condi- 
tions for achieving a phantom phase, listed above. Hence, 



we have verified that the late time attractor solution for 
V{4>) = is stable once we have [i such that we get a 
phantom phase. 

Now for generality we will repeat the same analysis for 
a negative potential, 



V((f>) = -Vbexp2a0. (A21) 
In this case the au- 



Hcrc wc will define y — — ' 
tonomous system will take the following form 



x' = -3.T + V6[m(1 + > 
+§£ [7 (l + x 2 + 
y' = V&axy + \y [7(1 



1 = 



— y 2 — x 2 



2 + y 2 )-ay 2 } + 
2 ) - 2x 2 ] 
x 2 + y 2 ) - 2x 2 ] 



(A22) 



There will be only two physical critical points in this case: 



r 2 VEfj, 

\ 3 -2+7 

1 1 V67 
\ 4 — a+n 



0) I.B 

1 \/2\/8^«-8^ 2 +67-37 2 , 
4 — a+fi > 



Just to clear any possible confusion please note that the 
labels I.B and II. B here are not related to the labeling 
of the different regions in Fig. 1. Notice that the first 
critical point in this case is identical to the third critical 
point for the case of positive potential, (see equation be- 
low (|A5[) ) The existence condition for the second critical 
point reads /i(/i — a) < §7(2 — 7). Fractional density 
values are 



Sid 



I.B 



, 3 !r 2 tS , "'" (A23) 

4 (-a+fj.) 2 



From here wc get an additional existence constraint on 
the second critical point, namely: 



a((j,-a) < -7 
The adiabatic constant for 6 is 



2a; 2 



14> 



37^ 



-4fia+4fj, 2 — 37 



I. B 

II. B 



(A24) 



(A25) 



Notice the difference between the last value in (|A25[) 
and the fourth in (|A9j) . Here we can set [i — > and 
recover 7 as the adiabatic index. Before the naive limit 
was —7 but one was not actually allowed to take that 
limit due to (|A7j) . Now the inequality has been reversed 
so the second critical point in the case of the negative 
potential could be used for a scaling solution in a phase 
where the hidden sector decouples from the ghost field. 

As before let us look at the effective equation of state 
parameter: 



Wtot 



1 3u h {u h ~l)+8p. 2 
3 

a— /x 



LB 
I I.B 



(A26) 
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Expanding around the critical points for the first solu- 
tion we get the following conditions for stability: 



> 



-8fi 2 - 3 7 2 + 12 7 - 12 



° > "2 



-2 + 7 
1 (-37 2 + 8pa - 8/i 2 + 67) 



-2 + 7 



(A27) 



The first inequality is trivially satisfied and the second 
implies: 



7 (2- 



7) < ^/^(m 



a) 



(A28) 



This could be used as a constraint on the steepness of 
the potential in order to preserve the tracking behavior 
in the phantom phase. For the second critical point we 
have the following eigenvalues: 



ei (2 ) 



1 1 



3aj — 6a + 6/i + Bj 



4 a — /z 

e 2( 2) = ^-^(3o7-6a + 6fi-flJ) (A29) 

Here B2 is a combination of the parameters, the exact 
form of which we will not be using. A necessary condition 
for the stability of this critical point is: 



a(l - Lo h ) < 2/j, 
The tracking ratios can be easily computed: 



(A30) 



(-2+ 7 )2 1+ 



(-2 + 7H J 

— 4a' 2 +4afj. — 3f 



LB 



II. B 



In this appendix we have found the conditions for scal- 
ing behavior in the phantom phase, for the cases of pos- 
itive, zero, and negative exponential potentials. An ex- 
ponential coupling between a ghost like field and some 
hidden sector fluid is used to generate a phantom phase, 
where the effective equation of state is less than negative 
one. We found the stable attractors and conditions nec- 
essary to obtain them. As a side result we notice that for 
negative potentials one could have a scaling stable late 
time attractor even if the ghost field is decoupled from 
the hidden sector. In that case the ghost will just track 
the hidden sector component. 



APPENDIX B: 



ESTIMATING THE MINIMUM 
VALUE OF M 



during the reheating phase (see Fig.l regions LA and 
LB). Once we enter the transition phase (region LA), if p, 
is large enough, most of the energy density will be stored 
in radiation, followed by the phantom field, and the least 
amount is contained by the hidden matter sector. To a 
good approximation we will use 'radiation domination' 
in what follows, since we are only interested in order of 
magnitude estimates. 

First we will look at the case when the hidden matter 
is non-relativistic when this transition occurs, i.e. lu = 0. 
Setting p r to zero, as appropriate for the transition to 
radiation, we will have the following equations for the 
evolution of the system: 



p h + 3Hp h = -p h T 
4> + 3H<j) = 



/3 7 + 4:Hpj = 

along with the Hubble equation 

1 







II 1 



3M| 



{Pi) 



(B2) 
(B3) 
(B4) 



(B5) 



We have neglected the annihilation term in equation {B4J) 
because p^T <C AHp 1 will be satisfied very quickly after 
the transition from the deSitter phase is started. The 
solutions for the energy densities are simple: 



Ph(t) 
p 7 (i) 



1 



t*Ci 
At 2 



2tm- 



(B6) 
(B7) 



where C\ is an integration constant we still have to fix. 
We have used the fact that during radiation domination 
H i in order to get the coefficient for the radiation 
energy density. Since it takes a longer time for the an- 
nihilation term to be sub-dominant with respect to the 
Hubble dampening term in the hidden matter continuity 
equation, we will have a short period where the energy 
density in hidden matter is decreasing faster than a~ 3 . 
This yields a minimum value for the parameter Sim- 
ply setting <%'(t) = leads to T(t c ) ~ H(t c ). In order to 
fix the integration constant C\ , we need to know at what 
time we start the transition from deSitter to radiation 
dominated phase. Using the initial condition 11 in (|3"5|) 
along with equation (|B7|) , the transition time is found to 
be: 



tlA 



Mlfc 
rrfirf 



(B8) 



Here we will estimate the minimum value of the ratio 
between the hidden sector energy density and the radia- 
tion energy density 

&=?±. (Bl) 
p 1 



Here we consider the transition from dS phase to the radiation 
dominated phase due to the sudden drop in the coupling between 
the hidden sector fields and the ghost field. The asymptotic value 
(valid for the dS phase) for the energy density in equation (1356 
(the right hand side of the equation) here becomes an initial 
condition for the transition phase, Region IA. 



15 



Furthermore, C\ is obtained by requiring that the hid- 
den matter density ph value in the deSitter phase matches 
the value estimated at the transition time tiA, obtained 
from (|B6[) . The initial conditions (|33[) together with (|3"2")) 
allows us to evaluate the initial value of hidden matter 
density at the beginning of the transition phase. In the 
limit of p p 3> 1 we get the following simplified form for 
the integration constant: 



Ci = 



V3(3M 2 + 32 v / 3/4m 3 ii j4 ) 



(B9) 



In order to get the minimum value for we go back to 
the condition T = Ph ^' ~ H(t c ), where t c represents 
the time at which this minimum is attained. This can be 
rewritten as: 



Ph 



Pi 



leading to 



V3M P 



2m 3 t c 
3M 2 



(BIO) 



(Bll) 



Also, from (|B10[) we can solve for t c using (|B6[) and (|B7[) : 
Plugging into (|B1 1|) we get the simplified 



tc = 

form: 



16 



32 



3M 2 m 3 C 2 



(B12) 



Using (|B8[) and (|B9|) in the above equation we obtain the 
final result, expressed only in terms of p p : 



<W ■ — 



1G 



V3/x|(l + 4VM^) 2 



(B13) 



Let us now turn our attention to the case where the 
hidden matter becomes relativistic at the energy scales 
where the transition between the deSitter and the radia- 
tion phases occurs. At the beginning there will still be a 
regime where the annihilation is effective, thus lowering 
the value of the ratio between the hidden matter energy 
density and the radiation energy density. If p r = 0, in- 
stead of a minimum we will now have a constant asymp- 
totic value towards which this ratio will tend. This is 
due to the fact that once the conversion is no longer effi- 
cient, both hidden matter and radiation energy densities 
will scale as a -4 . For completeness and generality we will 
study the case where p r ^ 0. As wc shall shortly see here 
a minimum develops, just as we have seen in the case of 
non-relativistic hidden matter. The approximative equa- 
tions we need to solve are: 



p h + AHp h 
4> + 3H<p 



-p h T + 2p h (p 



Pr 



(B14) 

(B15) 
(B16) 



along with the Hubble equation 

1 



H< = 



3M; 



-p~ 



(B17) 



Notice that we have neglected the ph terms in the second 
and third equations since we are interested in the phase 
of rapid conversion of hidden sector particles to radia- 
tion. Therefore, in this regime the energy density stored 
in the hidden sector will decay much faster than a~ 4 . Es- 
sentially the reason for this being that p has transitioned 
from large to small values, i.e p p 3> p r . The system 
admits the following solutions: 



H(t) 

3A.„ 
Py(t) = - 

m = 



i 

~2t 

At 2 



8m 3 A (j> 2 p r 2 



Ph (*) 



Where D(t) is: 



(B18) 
(B19) 
(B20) 

(B21) 

,3 a 2,, 2 



D(t) = Ar^MpA^pr + rM/ + 8t 2 e M ^C h m 6 A/p 

In order to obtain the integration constants Ch and A^ 
we match the energy densities above, evaluated at the 
transition time t\A with the corresponding values from 
the deSitter phase, given in f|32|) to l|35p . As before, t\A 
is given by (|B8j) . For A^ this leads to: 



2 p p M p 2 ^ 4 e /4 

3 ^5/4 m 3/2 



(B22) 



The expression for Ch for general p p turns out to be 
messy, but it can be simplified in the limit p p S> 1 to: 

1 V3p p {-Ap r 2 +Ap p p r +p p )c- 4 ^ 

Ch = -? 771 — 2 t^ 3 ' 

3 M p p r z 

Above we have used £ ~ anc j ^ ^ 4\/3/Up, expres- 
sions valid in the p p 3> 1 limit. 

Since we have p r ^ it we expect a minimum of the 
ratio $% to develop. Re-expressing St = by use of the 
definition of the ratio between the hidden sector and 
the radiation energy densities, and of the equations (|B14|) 
and (|B16|) we get, 



r (ph + p r ) = 2pr4> 



Pr 



(B24) 



Defining t c the time at which the above equation holds, 
we find that 



2m 3 (j)(t c )p r 8 m 3 'p r A^y/t 
M p p 1 (t c ) 3 



M p 3 



(B25) 
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From (|B24[) we can also obtain t c by using the solutions 
we have for the energy densities in (|B19|) and (|B21|1 . In 
order to get a closed form we will need to do some ap- 
proximation of the exponent in the denominator of (|B21[) . 
Since t > t\A we can show using t\A from (jB8[) and 
from (|B22|) that 3^7= < 4^ r . Since /j r « 1 wc will trun- 
cate the expansion of the exponential at terms of 0([i 2 ), 
leading to: 



C h 2 m 6 A^ 6 fx, 6 



t r = 1024 ■ 



(B26) 



M 2 (M p 2 + 8 CW VMr 2 ) 
Going back to (|B25j) we get, after some algebra, the fol- 



lowing form for & m 



a*. -A^! e " 4 ^ 

t y/ mm r— o c 

V3 /V 



- (i p (flp + 4Hr(flp - Mr)) 1 



We can double expand this expression in ^ r and [i p 1 and 
keep only the leading terms: 



1 V3(3 + 8^r) 



+ 0(ii 2 r ,^ 3 ) 
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In the limit /i r — > we get the result in eq.([4T 



(B27) 
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